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A  laser  beam  propagating  in  the  atmosphere  undergoes  scattering  and  absorption  by 
atmospheric  molecules  and  suspended  particles.  A  portion  of  the  scattered  radiation  is 
backscattered  and  propagates  toward  the  location  of  the  laser  source.  There  is  concern 
that  the  backscattered  radiation  produced  by  a  high-energy  laser  may  pose  ocular  or  skin 
hazards  to  personnel  in  the  vicinity  of  the  laser  source.  To  assess  these  hazards,  equations 
are  derived  for  the  backscattered  irradiance  from  both  continuous-wave  and  pulsed 
lasers.  The  derivations  are  based  on  the  definition  of  the  atmospheric  backscatter 
coefficient  and  the  principles  of  beam  propagation  theory.  The  expression  for  the 
backscatter  irradiance  of  a  pulsed  laser  is  shown  to  reduce  to  the  LIDAR  ( Light  Detection 
And  Ranging)  equation  under  conditions  that  characterize  typical  remote  sensing 
scenarios.  The  application  of  the  derived  equations  to  backscatter  hazard  assessments 
using  the  applicable  ANSI  standards  for  safe  use  of  lasers  is  discussed. 
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Nomenclature 

limiting  aperture 

area  element  in  the  x-y  plane  with  its  center  at  the  origin 
volume  element  in  the  space  illuminated  by  the  beam 
differential  light  scattering  cross  section 

solid  angle  subtended  by  d A  at  the  center  of  the  volume  element  dV 
continuous-wave  (CW)  beam  irradiance 
pulse  irradiance 

CW  beam  backscattered  irradiance  at  the  center  of  the  laser  aperture 
CW  beam  backscattered  irradiance  on  the  x-y  plane  at  distance  p  from  the 
center  of  the  laser  aperture 

backscattered  pulse  irradiance  at  the  center  of  the  laser  aperture  at  time  t 
backscattered  pulse  irradiance  on  the  x-y  plane  at  time  t  and  distance  p  from 
the  center  of  the  laser  aperture 
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EsX  irradiance  of  a  LIDAR  return  pulse 
/  beam  focus  distance 

hp(z)  altitude  at  the  location  on  the  slant  propagation  path  that  is  at  a  distance  z  from 
the  aperture  center 
M2  beam  quality 

Nk  number  of  scatterers  of  type  k  per  unit  volume 

Q  pulse  energy 

<2j,  p  effective  energy  of  backscattered  pulse 

/8(z)  backscatter  coefficient  at  distance  z  from  the  aperture  center  along  the  propa¬ 
gation  path 

9b  beam  1/e  half-angle  divergence 

X  laser  wavelength 

/r  atmospheric  extinction  coefficient 

|  normalized  scalar  field  in  the  region  z  >  0 

p  distance  from  the  origin  on  the  x-y  plane 

r  pulse  duration 

ra(z)  atmospheric  attenuation  from  the  aperture  center  to  location  z  along  the  prop¬ 
agation  path 

<t>  CW  beam  power 

d>j,cw  effective  power  of  the  backscattered  radiation  produced  by  a  CW  beam 
backscattered  pulse  power  passing  through  the  limiting  aperture  at  time  t 
cd0  beam  Me2  radius  on  the  laser  aperture 


1.  Introduction 

Laser  radiation  propagating  in  the  atmosphere  is  scattered  and  absorbed  by  molecules 
and  suspended  solid  particles  and  liquid  droplets.  For  high-energy  lasers,  the  scattered 
radiation  may  pose  ocular  or  skin  hazards  that  must  be  quantified  and  evaluated  prior  to 
field  tests  or  operational  deployment.  Of  particular  concern  is  the  backscattered  radiation, 
i.e.,  the  radiation  scattered  in  the  direction  opposite  to  the  beam  propagation  direction, 
which  may  pose  a  hazard  to  personnel  in  the  vicinity  of  the  laser  source. 

Backscattered  laser  radiation  is  commonly  used  to  measure  atmospheric  properties  with 
a  technique  known  as  LIDAR  (Light  Detection  And  Ranging).  In  a  typical  LIDAR  applica¬ 
tion,  laser  pulses  are  transmitted  vertically  into  the  atmosphere  and  the  backscattered  signal 
is  processed  using  the  range  resolution  capabilities  to  estimate  atmospheric  properties  as  a 
function  of  altitude.7  LIDAR  systems  have  been  used  to  measure  such  atmospheric  proper¬ 
ties  as  temperature,  pressure,  humidity,  wind  speed  and  direction,  and  aerosol  (suspended 
particles  and  droplets)  and  trace  gas  concentrations.3,4 

In  LIDAR  systems,  the  information  of  interest  in  the  backscattered  signal  originates 
at  long  distances  from  the  laser  source,  typically  at  least  a  kilometer.  This  information 
therefore  arrives  at  the  receiver  some  time  after  the  outgoing  pulse  left  the  laser  aperture. 
However,  because  backscattered  laser  radiation  is  produced  as  soon  as  the  beam  starts 
to  exit  the  laser  aperture,  hazard  assessments  must  quantify  the  backscattered  irradiance 
starting  at  this  time. 

This  paper  provides  a  quantitative  framework  for  hazard  assessments  of  the  backscattered 
laser  radiation  produced  by  both  continuous-wave  (CW)  and  pulsed  beams.  Models  for 
CW  and  pulsed  beams  are  derived  and  then  applied  to  obtain  approximate  expressions 
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that  quantify  the  backscattered  irradiance  at  the  plane  of  the  laser  aperture.  The  derived 
expressions  for  the  backscattered  irradiance  of  a  pulsed  beam  are  valid  at  all  times  after  the 
pulse  emerges  from  the  laser  aperture.  Unlike  the  backscattered  irradiance  produced  by  a 
CW  beam,  which  does  not  vary  in  time,  the  backscattered  irradiance  produced  by  a  pulsed 
beam  is  shown  to  exhibit  pulse-like  characteristics. 

A  hazard  assessment  of  the  backscattered  radiation  from  a  CW  or  a  pulsed  laser  beam 
requires  the  calculation  of  several  parameters  to  implement  the  hazard  evaluation  procedures 
specified  in  the  ANSI  standard  for  safe  use  of  lasers.1  We  discuss  how  the  expressions  for 
the  backscattered  irradiance  derived  in  this  paper  are  applied  to  calculate  the  required 
parameters. 


2.  Backscatter  Coefficient 

J  ' 

The  backscatter  coefficient  is  a  measurable  quantity  that  determines  the  backscattered 
irradiance.  In  this  section,  the  backscatter  coefficient  is  used  to  derive  the  differential 
backscattered  irradiance  at  the  center  of  the  laser  aperture.  This  irradiance  is  assumed 
to  originate  in  a  small-volume  element  illuminated  by  the  beam.  The  total  backscattered 
irradiance  at  the  aperture  center  is  then  found  by  integrating  the  backscattered  irradiance 
contributions  of  individual  volume  elements.  In  Sec.  3  we  present  beam  models  that, 
together  with  the  expression  for  the  differential  irradiance  derived  in  this  section,  are  used 
to  derive  the  total  backscattered  irradiance  for  both  CW  and  pulsed  beams. 

The  differential  light  scattering  cross  section,  dcr(0)/d£2,  of  a  single  molecule  or  particle 
has  units  of  square  meters  per  steradian  and  relates  the  intensity  of  the  light  scattered  at 
angle  6, 1(6)  in  units  of  watts  per  steradian  to  the  incident  irradiance  E: 


1(6)  =  E  ■ 


d  o(6) 
df2  ' 


(1) 


The  scatter  angle  6  is  measured  counterclockwise  from  the  propagation  direction  of  the 
incident  beam  so  that  at  6  =  it  (6  =  0),  the  reflected  intensity  is  in  the  direction  opposite 
(parallel)  to  the  beam  propagation  direction.  In  general,  the  differential  scattering  cross 
section  of  a  single  molecule  or  particle  depends  on  the  laser  wavelength. 

The  backscatter  coefficient,  f),  is  defined  by  the  expression 


#=£ 

k 


dcr^TT) 

dfi 


■Nk, 


(2) 


where  Nk  is  the  number  of  scatterers  of  type  k  per  unit  volume.  The  backscatter  coefficient 
consists  of  contributions  from  all  molecular  and  aerosol  species  that  are  present  in  the  atmo¬ 
spheric  volume  illuminated  by  the  beam.  It  follows  from  Eq.  (2)  that  /3  has  units  of  meters 
per  steradian  and,  in  general,  depends  on  altitude,  the  prevailing  atmospheric  conditions, 
and  the  wavelength  of  the  incident  beam  irradiance.  We  note  that  the  molecular  contribution 
to  /3  is  due  primarily  to  the  nitrogen  and  oxygen  that  make  up  99%  of  the  atmosphere’s 
molecular  constituents.  The  differential  light  scattering  cross  section  of  molecules  is  given 
by  the  Rayleigh  formula.5 

Figure  1  shows  the  coordinate  system,  the  irradiance  distribution  of  a  hypothetical 
focused  CW  beam,  and  the  differential  volume  and  area  elements  used  in  the  derivation  of 
the  differential  backscattered  irradiance.  The  laser  aperture  lies  in  the  x-y  plane,  and  the 
beam  is  assumed  to  propagate  in  the  positive  z  direction  along  an  arbitrary  slant  path  through 
the  atmosphere.  This  path  begins  at  the  coordinate  system  origin  located  at  the  center  of 


Journal  of  Directed  Energy ,  3,  Summer  2010 


292 


MEGALOUDIS,  EARLY,  AND  KENNEDY 


Fig.  1.  Focused  laser  beam  and  the  differential  area  and  volume  elements  used  to  derive 
the  backscattered  irradiance. 


the  laser  aperture,  so  that  the  z  axis  represents  the  distance  from  the  origin  along  the  path. 
Because  both  the  backscatter  and  atmospheric  extinction  coefficients  vary  with  altitude, 
they  will  also  vary  over  the  propagation  path  unless  it  is  at  constant  altitude.  To  account  for 
this  variation  with  path  distance,  we  introduce  the  path-dependent  backscatter  coefficient, 
/6(z),  and  atmospheric  attenuation,  ra(z).  Both  of  these  quantities  may  be  obtained  from 
backscatter  and  extinction  measurements  or  model  predictions  at  different  altitudes.  For 
example,  let  represent  the  backscatter  coefficient  as  a  function  of  altitude  h.  Also, 
let  the  function  hp(z)  represent  the  altitude  at  the  location  on  the  slant  path  that  is  a 
distance  z  from  the  origin.  Then,  the  path-dependent  backscatter  coefficient,  /3(z),  is  given 

by  A nlhpiz)). 

Consider  the  differential  volume  element  dV  shown  in  Fig.  1,  defined  as  dV  =  dxdydz 
and  centered  at  location  (x,  y,  z).  Using  the  definition  of  /3  and  the  beam  irradiance  incident 
on  dV ,  we  find  that  the  differential  backscattered  intensity,  dls,  at  the  center  of  dV  is  given 
by 


d/s  =  /J(z)  ■  E{x,  y,  z)  ■  ra(z)  •  dx  dy  dz, 


(3) 


where  E  is  the  unattenuated  beam  irradiance  at  dV  and  za(z)  accounts  for  atmospheric 
attenuation  of  the  beam  over  the  propagation  path  from  the  aperture  to  dV .  The  solid  angle 
d£2,  subtended  by  dA  at  the  center  of  dV,  is 


d£2  = 


Z  •  dA 

(x2  -I-  y2  +  z2)1-5 ' 


(4) 


Therefore,  the  backscattered  power,  dP,  passing  through  the  area  dA  centered  at  the  origin 
as  shown  in  Fig.  1  is  approximated  by  multiplying  dls  by  dQ  so  that  dP  is  given  by 


dP  =  ■  E(x,  y,  z)  • 


Z  •  dA 


(x2  +  y2  +  z2) 


2\1.5 


r^(z)  •  dx  dy  dz. 


(5) 
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In  this  expression,  we  have  included  an  additional  factor  of  rfl(z)  to  account  for  atmospheric 
attenuation  over  the  propagation  path  from  dV  back  to  the  x-y  plane.  Dividing  both  sides 
of  Eq.  (5)  by  dA,  we  find  that  the  differential  backscattered  irradiance  at  the  origin  is 
approximately  given  by 

d  Es  =  fi(z)  ■  E(x,  y,  z)  ■  2  \  2  ,,  •  r2(z)  •  dx  dy  dz.  (6) 

(x2  +  y2  +  z2)1-5 

From  Fig.  1,  one  can  see  that  for  volume  elements  near  the  x-y  plane  and  relatively  far 
from  the  z  axis,  the  scatter  angle  6  to  the  origin  will  not  be  n  radians  as  assumed  in  the 
derivation  of  Eq.  (6).  However,  for  the  great  majority  of  locations  in  the  beam,  the  angle 
to  the  origin  will  be  close  to  n  radians  because  the  extents  of  the  beam  along  the  x  and  y 
directions  are  small  compared  to  its  extent  along  the  z  direction.  Moreover,  for  locations 
near  the  aperture  plane  and  away  from  the  z  axis,  the  subtended  solid  angle  is  small  so  that 
their  contribution  to  the  irradiance  at  the  origin  will  also  be  small. 

Equation  (6)  shows  that  in  order  to  assess  the  backscattered  hazard  associated  with  a 
specific  high-energy  laser  test  or  deployment,  one  needs  to  know,  or  at  least  be  able  to 
estimate  an  upper  bound  for,  the  backscatter  coefficient  at  the  laser  wavelength  and  the  test 
or  deployment  location.  This  is  complicated  by  the  variability  of  aerosol  species  and  their 
concentrations.  Aerosols  also  vary  in  shape,  with  the  majority  having  radii  between  10-3 
and  10  fxm.  Examples  of  aerosols  include  particles  of  dust,  combustion  by-products,  fog, 
and  possibly  salt  when  droplets  of  sea  spray  evaporate.  The  scattering  of  light  by  spherical 
particles  with  arbitrary  refractive  index  is  described  by  the  theory  developed  by  Mie.2 
However,  because  aerosols  are  rarely  spherical,  the  Mie  theory  is  not  always  accurate. 

Despite  these  challenges,  models  of  aerosol  light  scattering  properties  have  been  devel¬ 
oped  based  on  measurements  as  well  as  theory.  These  models  are  representative  of  rural, 
urban,  maritime,  and  tropospheric  atmospheric  conditions  and  are  documented  in  Ref.  6. 
For  each  model,  this  document  gives  coefficients  of  extinction,  scattering,  and  the  angular 
scattering  distribution  as  functions  of  wavelength  and  relative  humidity.  In  the  absence  of 
site-specific  backscatter  measurements,  the  data  may  be  used  to  estimate  an  upper  bound 
for  the  backscatter  coefficient  at  the  location  of  interest. 

3.  Beam  Irradiance  Models 

According  to  Eq.  (6),  calculation  of  the  backscattered  irradiance  requires  that  we  know  the 
beam  irradiance  distribution  in  the  scattering  volume.  In  this  section,  we  present  parametric 
beam  irradiance  models  that  can  be  used  to  predict  the  backscattered  irradiance  from  a  CW 
or  pulsed  laser  beam  based  on  knowledge  of  beam  parameters. 

The  unattenuated  irradiance  of  a  CW  beam  propagating  parallel  to  the  z  axis,  few.  is 
modeled  as 

£cw(*,  y.z)  =  $  •  lt(*,;y,z)l2,  (7) 

where  <t>  is  the  beam  power  and  £(x,  y,  z)  is  the  field  in  the  space  to  the  right  of  the  x-y 
plane  as  shown  in  Fig.  1. 

The  pulsed  beam  model  assumes  that  all  pulses  in  the  beam  are  identical  and  that  they 
propagate  along  the  positive  z  axis.  With  these  assumptions  the  unattenuated  irradiance, 
Ep,  of  a  single  propagating  beam  pulse  is  modeled  as 

Ep(x,y,z,t)=  ®  p  (t  -  •!£(■*, y,z)\2,  (  (8) 
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where  Q  is  the  pulse  energy  and  r  is  the  pulse  duration.  The  dimensionless  function  p(t) 
is  defined  as 


pit )  =1  for  0  <  t  <  r, 

pit)  =  0  otherwise.  (9) 


Therefore,  at  a  given  time  t,  the  leading  and  trailing  edges  of  the  pulse  are  at  locations 
z  =  ct  and  c(t  —  r),  respectively.  It  follows  that  at  t  =  0  the  leading  edge  of  the  pulse  is  at 
z  =  0,  i.e.,  in  the  plane  of  the  exit  aperture. 

The  field  £(x,  y,  z)  is  normalized  so  that  for  any  value  of  z,  its  integral  over  the  x-y 
plane  is  unity: 

r  f°°  \Hx,y,z)\2dxdy=  1.  (10) 

J- 00  J -00 


This  condition  ensures  that  the  integral  of  Eq w  over  any  plane  perpendicular  to  the  prop¬ 
agation  direction  equals  the  beam  power.  Equation  (10)  may  also  used  to  show  that  the 
integral  of  the  single  pulse  irradiance  given  by  Eq.  (8)  over  the  pulse  duration  and  the  plane 
perpendicular  to  the  propagation  direction  equals  the  pulse  energy,  i.e., 

p(z/c)+T  poo  poo 

/  At  /  /  EP(x,y,z,t)dxdy  =  Q.  (11) 

J  z/c  J —00  J —00 

In  general,  the  beam  spatial  irradiance  distribution  will  depend  on  the  output  from  the 
laser  and  the  optics  used  to  transmit  the  beam.  We  use  a  beam  irradiance  model  that  has  a 
Gaussian  distribution,  accounts  for  the  effect  of  higher  order  modes  that  may  be  present, 
and  describes  unfocused  as  well  as  focused  beams.  In  this  model,  |£(jc,  y,  z) \2  is  given  by 


!£(*,?,  z)l2 


exp  [-(x2  +  y2)/co2iz)] 
7t  ■  <W2(z) 


where  o>(z)  is  the  He  beam  radius  and  o>2(z)  is  given  by 


(12) 


( 


X 

n  ■  co0 


+  coi 


(13) 


Here,  a>0  is  the  1/e2  beam  radius  at  the  aperture.  This  model  satisfies  Eq.  (10)  and  was 
derived  with  the  assumption  that  the  beam  is  not  truncated  by  the  aperture.  In  Eq.  (13),  X 
is  the  laser  wavelength,  /  is  the  focus  distance,  and  M2  is  the  beam  quality.  If  the  focus 
distance  /  goes  to  infinity  in  Eq.  (13),  the  model  describes  an  unfocused  beam  with  its 
waist  at  the  aperture  and  radius  co0.  In  this  case,  the  beam  diverges  with  1/e  half-angle  9b 
given  by 


n  ■  (lig 


(14) 


which  is  seen  to  be  M 2  times  the  1/e  half-angle  divergence  of  the  lowest-order  Gaussian 
mode.  In  the  literature,  the  divergence  of  this  mode  is  often  referred  to  as  the  “diffraction 
limit.”  For  this  reason,  some  authors  refer  to  M 2  as  the  “times-diffraction-limit”  parameter. 
The  inclusion  of  beam  quality  in  the  model  makes  it  possible  to  approximately  account  for 
the  presence  of  higher  order  modes  that  have  a  higher  beam  divergence  than  the  lowest-order 
Gaussian  mode. 
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Note  that  the  expressions  for  the  backscattered  irradiance  derived  in  the  next  section 
are  also  valid  for  other  irradiance  models  that  are  rotationally  symmetric  about  the  beam 
propagation  axis.  These  models  include,  for  example,  top-hat  beams  and  centrally  obscured 
beams  transmitted  by  on-axis  telescopes.  To  derive  the  backscattered  irradiance  for  any  of 
these  models,  one  has  only  to  replace  our  model  irradiance  distribution  [Eqs.  (12)  and  (13)] 
with  that  of  the  new  model  in  the  integrals  for  the  backscattered  irradiance.  The  model 
irradiance  distribution  must  be  normalized  in  accordance  with  Eq.  (10). 


4.  Backscattered  Irradiance 

4.1.  CW  beam 


The  total  backscattered  irradiance  at  the  origin  is  found  by  summing  the  backscatter 
contributions  from  individual  volume  elements  in  the  space  occupied  by  the  beam,  or 
equivalently  by  integrating  Eq.  (6).  Substituting  the  CW  model  irradiance  given  by  Eq.  (7) 
in  Eq.  (6)  and  integrating  the  result  gives  the  total  backscattered  irradiance  at  the  origin, 
Es.cw'- 


'S.CW 


rOO  noo  rot 

=  *■/  /  / 

JO  J  —  oo  J  —  o 


lg(x,  y.  z)l2  •  P(z)  ■  r2(z)  ■  Z 
(. x 2  +  y2  +  z2)1-5 


dx  dy  dz. 


(15) 


This  integral  may  be  simplified  by  exploiting  the  rotational  symmetry  of  || (x,  y,  z)|2  about 
the  z  axis.  If  we  define  | £(r,  z)|2  as 


2  _  exp  [— r2/&>2(z)] 


2/  \ 

7t  ■  a/(z) 

and  change  the  x  and  y  variables  of  integration  by  the  transformation 


(16) 


x  =  r  ■  cos  <p, 
y  =  r  ■  siny, 


(17) 


we  find  that  Eq.  (15)  may  be  written  as 


ij.CW 


=  2  •  n  ■  4> 


r  OO  /» o 

Jo  Jo 


I f(r,  z)f  •  &(z)  •  ra2(z)  •  Z 
(ri  +  z2)1-5 


r  ■  dr  dz. 


(18) 


The  numerical  evaluation  of  Eq.  (18)  is  faster  than  that  of  Eq.  (15). 

Note  that  if  the  propagation  path  is  one  of  constant  altitude,  /i0,  then  fi{z)  in  the  integrand 
of  Eq.  (18)  is  a  constant  equal  to  the  atmospheric  backscatter  coefficient  at  altitude  ho  and 
may  therefore  be  taken  out  of  the  integral.  Also,  for  this  propagation  path  the  two-way 
atmospheric  attenuation  in  the  integrand  is  given  by  exp[— 2  •  fx(ho)  ■  z],  where  fi(ho)  is  the 
atmospheric  extinction  coefficient  at  altitude  ho. 


4.2.  Pulsed  beam 


Consider  the  differential  volume  element  dV  shown  in  Fig.  1  and  a  single  pulse  that 
propagates  in  the  positive  z  direction.  The  power  backscattered  at  dV  and  time  t  into  the 
solid  angle  d£2  is  denoted  by  dJfoy(f)  and  is  given  by 


dPdv(f)  =  Kz)  ■ 


EP(x,y,z,t)  ■  dA  -  z  •  t a(z) 
(x2  +  y2  +  z2)1-5 


dx  dy  dz, 


(19) 
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where  EP(x,  y,  z,  t)  is  the  pulse  irradiance  given  by  Eq.  (8).  The  propagation  from  dV 
to  the  aperture  plane  introduces  an  additional  time  delay,  so  the  backscattered  power  that 
passes  through  dA  is 


dPdx(f)  =  ftz)  • 


Ep[x,  y,z,t  —  (z/c)]  • dA  •  z  •  r2(z) 
(x2  +y2  +  z2)1-5 


dx  dy  dz, 


(20) 


where  we  have  included  an  additional  factor  of  ra(z)  to  account  for  the  attenuation  over  the 
path  back  to  the  x-y  plane.  If  we  make  use  of  Eq.  (8)  in  Eq.  (20)  and  divide  by  dA,  we  find 
that  the  differential  backscattered  irradiance  at  the  origin,  dEsP(t),  is  given  by 


d E,At)  =  P(z)  •  - 
x 


I5(*.  y,z)l2  •  z  ■  x 2(z) 
( x 2  +  y2  +  z2)1-5 


•  dx  dy  dz. 


(21) 


At  a  fixed  time  f,  the  factor  p[t  —  (2 z/c)]  is  nonzero  for  values  of  z  such  that  0  < 
[r  —  (2z/c)]  <  r  or  equivalently  if  z  lies  between  [c(r  —  r)]/2  and  cf/2.  In  other  words, 
the  irradiance  at  the  origin  at  an  arbitrary  time  t,  dEs  P(t),  is  determined  when  the  leading 
and  trailing  edge  of  the  pulse  is  located  at  ct/2  and  [c(?  —  r)]/2,  respectively.  Therefore, 
the  total  backscattered  irradiance  at  the  origin  at  time  t  is  obtained  by  integrating  Eq.  (21) 
to  yield 


EsAD 


Q  ret/ 2  roo  * 

X  J c(t— r)/2  J—  oo  J— oo 


I5(*.y.  z)l2  ■ /8(z)  •  r2(z)  ■  z 
(x2  +  y2  +  z2)1-5 


dx  dy  dz.  (22) 


Equation  (22)  is  valid  when  t  >  x.  However,  for  times  t  such  that  0  <  t  <  r,  the  pulse 
trailing  edge  location,  [c(r  —  r)]/2,  will  be  less  than  zero,  which  means  that  some  portion 
of  the  pulse  is  in  the  region  to  the  left  of  the  exit  aperture  plane.  This  portion  of  the  pulse 
does  not  contribute  to  the  backscattered  irradiance  at  the  origin,  so  that  for  times  t  such 
that  0  <  t  <  x,  the  lower  limit  in  the  integral  over  z  must  be  set  to  zero.  Therefore,  the 
backscattered  irradiance  at  the  origin  from  a  single  outgoing  pulse  is  given  for  all  times 
t  >  0  by 


Q  ret/2  poo  p 

EsAO  =  -■  /  / 

X  J  a  J  —oo  J  — 


l£(x,  y,  z)|2  •  P(z)  ■  x2(z)  ■  z 
(x2  +  y2  +  z2)1-5 


dx  dy  dz. 


(23) 


where 


a  =  0  for  0  <  t  <  r, 


a  — 


c(t  -  x) 


for  t  >  x. 


(24) 


The  integral  in  Eq.  (23)  may  be  simplified  by  making  use  of  the  rotational  symmetry  of 
||(x,  y,  z)|2  about  the  z  axis.  Changing  variables  of  integration  as  in  Eq.  (17),  we  find  that 


Es,p(t )  =  2  •  7T 


q  r"1  r  fr.  z)\2  ■  m  •  t2(Z)  •  z 
T'i  /„  - - rdrdz 


(25) 


We  note  that  the  ratio  of  the  pulse  energy  to  the  pulse  duration  is  called  peak  power. 
Equation  (25)  shows  that  the  backscattered  power  is  proportional  to  the  product  of  the 
backscatter  coefficient  and  the  peak  power.  Finally,  the  backscattered  pulse  radiant  exposure 
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Fig.  2.  Typical  backscattered  irradiance  waveform  on  the  aperture  plane  produced  by  a 
propagating  laser  pulse. 


at  the  origin,  HSiP,  is  given  by 

/•OO 

HS,P=  ES'P(t)-dt.  (26) 

Jo 

Figure  2  shows  the  typical  temporal  behavior  of  the  backscattered  irradiance  from  a  sin¬ 
gle  outgoing  pulse  as  calculated  with  Eq.  (25).  For  this  calculation,  we  assumed  a  horizontal 
propagation  path,  a  laser  wavelength  of  1 .064  /urn,  a  beam  radius  at  the  aperture  of  0.2  m,  a 
pulse  energy  and  duration  of  1  J  and  50  ns,  respectively,  and  a  beam  focus  distance  of 
1  km.  The  backscatter  and  extinction  coefficients  were  2  x  10-6  m-1-  sr-1  and  9  x 
10-6  m-1,  respectively.  For  other  input  parameter  values,  the  temporal  behavior  of  the 
backscattered  irradiance  is  qualitatively  similar  to  that  shown  in  Fig.  2. 

Figure  2  shows  that  for  times  between  zero  and  the  pulse  duration,  the  backscattered 
irradiance  rises  rapidly  from  zero  to  its  maximum  value.  After  reaching  its  maximum,  the 
backscattered  irradiance  decays,  so  that  after  a  time  of  500  ns  its  value  has  been  reduced 
by  a  factor  of  10-4  from  its  maximum  value.  From  Fig.  2  we  can  see  that  the  backscattered 
irradiance  resembles  a  pulse  waveform.  If  we  define  the  term  pulse  width  as  the  time  interval 
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over  which  a  pulse  is  reduced  to  half  its  maximum  value,  we  can  see  from  Fig.  2  that  the 
backscattered  pulse  duration  is  slightly  longer  than  the  width  of  the  pulse  in  the  beam. 

We  may  therefore  conclude  that  every  pulse  of  a  pulsed  beam  that  propagates  in  the 
atmosphere  will  generate  a  backscattered  pulse  that  propagates  back  toward  the  laser  source. 
Therefore,  an  observer  at  or  near  this  location  will  be  exposed  to  a  train  of  backscattered 
pulses.  In  Sec.  5,  we  discuss  how  the  results  of  this  section  can  be  applied  to  assess  the 
hazard  posed  by  the  backscattered  pulse  train. 


4.3.  Derivation  of  the  LIDAR  equation 


LIDAR  is  a  technique  in  which  a  pulsed  laser  is  used  to  make  remote,  range-resolved 
measurements  of  the  properties  of  the  atmosphere.  In  addition  to  a  transmitter  that  generates 
the  outgoing  beam,  a  LIDAR  system  employs  a  receiver  that  collects  the  return  pulses 
produced  by  transmitted  pulses  that  have  been  backscattered  by  atmospheric  aerosols  and 
molecules.  In  remote  sensing  applications,  the  collected  return  pulses  originate  from  regions 
that  are  at  large  distances  from  the  system.  Equivalently,  the  time  at  which  the  return  pulse 
arrives  at  the  receiver  is  very  much  longer  than  the  pulse  duration.  Under  these  conditions, 
the  irradiance  of  the  return  pulse,  Es  i(t),  is  given  by  the  LIDAR  equation 


£*,l(0  =  £(R) '  *2  '  \  •  '  exP 


-2-  f  n(z)dz  , 
Jo 


(27) 


where  R  =  ct/2  and  /i(z)  is  the  path-dependent  atmospheric  extinction  coefficient  defined 
in  a  manner  analogous  to  (${z).  The  last  factor  on  the  right  side  of  Eq.  (27)  is  the  square 
of  the  path-dependent  attenuation  at  z  =  R,  z2(R).  As  a  check  on  the  validity  of  Eq.  (23), 
we  now  show  how  the  LIDAR  equation  can  be  derived  from  Eq.  (23)  under  the  appropriate 
assumptions. 

At  times  t  that  are  much  longer  than  r,  the  pulse  will  have  traveled  a  distance  that  is  very 
much  greater  than  the  characteristic  radius  of  |£(jc,  y,  z) |2  in  the  x-y  plane.  The  quantity 
|f(x,  y>  z)\2,  and  hence  the  integrand  in  Eq.  (23),  will  be  small  for  x  and  y  outside  the 
region  bounded  by  this  radius.  Therefore,  the  sum  x2  +  y2  in  the  denominator  will  be  much 
smaller  than  z2  and  Eq.  (23)  may  be  approximated  by 


W)Sif  rr 

t  J  [c(r— r)]/2  J  — oo  J  — oo 


!£(*>  y,  z)\2  ■  fi(z)  ■  t2(z) 


dx  dy  dz. 


(28) 


Using  Eq.  (10)  we  can  perform  the  integration  over  x  and  y: 


EsAt ) 


q  rc,/2  m  ■  r 2(z)  d_ 

T  J[c(t- r)]/2  z2 


(29) 


The  change  in  fi  and  xa  over  the  limits  of  integration  is  negligible,  so  that  both  may  be 
taken  out  of  the  integral  and  evaluated  at  z  =  ct/2.  The  resulting  integral  is 


r'2  -  t 

[c(r-t)]/2  z2  (c/2)r  •  (t  -  r)’ 


(30) 


and  because  t  is  assumed  to  be  much  greater  than  r,  we  may  approximate  the  right-hand 
side  by  r/(c/2)  •  t2.  Substituting  these  results  in  Eq.  (29),  we  find  that 


EsAO  =  Pict/2)  Q  - 


(ct/2? 


■  exp 


fCI/2 

-2-  , 
Jo 


/x(z)  ■  dz 


(31) 
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which  is  Eq.  (27).  Thus,  we  have  shown  that  the  widely  used  LIDAR  equation  may  be 
derived  from  Eq.  (23)  under  assumptions  that  are  valid  for  remote  sensing  scenarios. 


5.  Application  to  Backscatter  Hazard  Assessment 

The  procedures  for  assessing  the  direct  beam  hazard  posed  by  CW  and  pulsed  laser 
beams  are  documented  in  the  ANSI  Z136.1-2007  standard.1  In  this  section,  we  discuss  how 
the  results  derived  earlier  in  this  paper  may  be  used  to  calculate  the  parameters  required  by 
these  procedures  in  order  to  evaluate  backscatter  hazards  from  CW  and  pulsed  beams. 


5.1.  CW  laser 


The  time-independent  backscattered  irradiance  produced  by  a  CW  beam  may  be  thought 
of  as  a  CW  beam  that  propagates  back  toward  the  x-y  plane,  potentially  creating  a  direct 
beam  hazard  for  an  observer  on  the  x-y  plane.  In  this  discussion  we  will  therefore  assume 
that  the  ANSI  methodology  for  evaluating  direct  CW  beam  hazards  can  be  applied  to  assess 
potential  hazards  associated  with  backscatter  from  a  CW  beam. 

For  a  CW  laser,  the  maximum  permissible  exposure  (MPE)  is  found  from  the  laser 
wavelength  and  the  total  exposure  duration.  In  addition  to  the  MPE,  the  procedure  specified 
in  the  ANSI  standard  also  requires  the  limiting  aperture  diameter  Df,  which  for  some  laser 
wavelengths  depends  on  the  exposure  time.  The  appropriate  limiting  aperture  diameter 
for  the  wavelength  and  exposure  time  of  interest  may  be  found  in  Table  8  of  the  ANSI 
standard.  For  example,  for  the  spectral  range  of  0.4 -1.4  /rm  and  for  exposure  durations 
between  10~13  and  3  x  104  s,  the  limiting  aperture  diameter  for  the  eye  is  7.0  mm.  The 
limiting  aperture  diameter  is  subsequently  used  to  calculate  the  effective  power,  which  is 
defined  to  be  the  power  transmitted  through  an  aperture  with  diameter  Df.  For  a  Gaussian 
beam  profile,  the  effective  power  is  calculated  from  a  closed-form  expression  that  depends 
on  Df  and  the  beam  diameter. 

An  analogous  closed-form  expression  does  not  exist  for  calculating  the  effective  power 
from  the  backscattered  irradiance.  To  derive  the  effective  power,  we  must  first  determine 
the  backscattered  irradiance  at  arbitrary  locations  on  the  x-y  plane.  This  is  accomplished 
by  exploiting  the  rotational  symmetry  of  our  beam  model  about  the  z  axis.  This  symmetry 
implies  that  the  backscattered  irradiance  will  be  constant  on  any  circle  in  the  x-y  plane 
whose  center  coincides  with  the  origin  of  the  coordinate  system  shown  in  Fig.  1.  In  other 
words,  the  backscattered  irradiance  will  vary  only  with  distance  from  the  origin  and  will 
be  independent  of  the  azimuth  angle  about  the  origin.  It  follows  that  the  variation  of  the 
irradiance  with  distance  from  the  origin  will  be  the  same  on  any  line  that  starts  at  the  origin. 
If  we  select  the  y  axis  as  this  line,  the  expression  for  the  backscattered  irradiance  as  a 
function  of  distance  from  the  origin,  p,  is  given  by 


r+oo  r+oo 

£j,cw(/0)  =  $>  ■  / 

Jo  Jo  Jo 


IC(a  z)I2  •  P(z)  •  ta2(z)  •  z 

[(r  •  cos  (p)2  +  (r  ■  sin  <p  —  p)2  +  z2]1-5 


r  ■  dtp  dr  dz. 

(32) 


We  should  note  that  a  plot  of  £s,cw  as  a  function  of  p,  which  can  be  generated  using 
Eq.  (32),  reveals  the  spatial  properties  of  the  backscattered  irradiance  on  the  x-y  plane. 
In  particular,  we  may  use  this  plot  to  visually  inspect  the  value,  width,  and  location  of 
the  maximum  backscattered  irradiance  on  the  x-y  plane.  A  representative  plot  of  the 
backscattered  irradiance  on  the  aperture  plane  as  a  function  of  distance  from  the  aperture 
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Fig.  3.  Example  backscattered  irradiance  on  the  aperture  plane  as  a  function  of  distance 
from  the  aperture  center. 


center  generated  using  Eq.  (32)  is  shown  in  Fig.  3.  For  this  calculation,  the  laser  wavelength 
was  1.064  gm ,  the  beam  radius  at  the  aperture  was  0.1  m,  the  beam  power  was  105  W, 
and  the  beam  focus  distance  was  1  km.  The  backscatter  and  extinction  coefficients  were 
2  x  10-6  m_1  sr-1  and  9  x  10-6  m-1,  respectively. 

As  shown  in  Fig.  3,  for  our  beam  irradiance  model  using  Eqs.  (12)  and  (13)  the  maximum 
irradiance  occurs  at  the  aperture  center.  If  the  maximum  backscattered  irradiance  occurs 
at  the  aperture  center,  the  effective  power  4>s, cw  is  found  by  integrating  E*,cw(p)  over  a 
circular  area  of  diameter  Df  centered  at  the  origin: 

r°r  n  _ 

4**,  cw  =  2  •  7r  •  /  £s,cw  (p)  •  P  d/x  (33) 

Jo 

A  hazardous  condition  exists  when  the  ratio  of  the  effective  power  to  the  product  of  the 
MPE  and  the  area  of  the  limiting  aperture  is  greater  than  unity.  In  the  event  that  Eq.  (32) 
has  its  maximum  at  a  point  on  the  aperture  plane  other  than  the  center,  then  the  effective 
power  must  be  calculated  by  integrating  ESi  cw(p)  over  a  circular  region  of  diameter  Df 
centered  at  the  point  of  maximum  irradiance. 
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If  a  hazardous  condition  exists  on  the  aperture  plane,  it  is  sometimes  useful  to  know  the 
distance  from  the  point  of  maximum  irradiance  to  a  location  at  which  the  irradiance  has  de¬ 
creased  to  a  nonhazardous  level.  To  find  this  distance,  the  calculation  of  the  effective  power 
for  different  circular  regions,  each  with  diameter  Df  but  centered  at  different  distances  from 
the  location  of  maximum  irradiance,  must  be  repeated.  The  effective  power  calculated  in 
each  case  will  decrease  with  increasing  distance  from  the  location  of  maximum  irradiance. 
The  distance  beyond  which  a  hazard  does  not  exist  will  be  that  for  which  the  corresponding 
effective  power  is  equal  to  the  product  of  the  MPE  and  the  area  of  the  limiting  aperture. 

Note  that  if  the  change  in  ESi  cw(p)  between  p  —  0  and  p  =  0.5  Df  is  small,  the  effective 
power  given  by  Eq.  (33)  can  be  approximated  by  the  product  of  the  limiting  aperture  area 
and  the  maximum  backscattered  irradiance,  £s,cw(0).  Therefore,  when  this  approximation 
is  valid,  the  ratio  of  the  effective  power  to  the  product  of  the  MPE  and  the  limiting  aperture 
area,  which  is  used  to  test  for  the  existence  of  a  hazardous  condition,  is  given  by  the  ratio 
of  the  maximum  backscattered  irradiance  to  the  MPE. 


5.2.  Pulsed  beam 


Based  on  the  discussion  of  Sec.  4.2,  each  pulse  that  propagates  in  the  atmosphere  will 
produce  a  backscattered  pulse  that  propagates  toward  the  location  of  the  laser  source. 
Therefore,  an  observer  at  or  near  this  location  will  be  exposed  to  a  train  of  backscattered 
pulses.  The  backscattered  pulse  train  is  analogous  to  that  produced  by  a  repetitively  pulsed 
laser.  In  this  discussion  we  will  therefore  assume  that  the  ANSI  methodology  for  assessing 
direct  beam  ocular  and  skin  hazards  from  repetitively  pulsed  lasers  can  be  applied  to  assess 
the  hazards  posed  by  the  backscattered  pulse  train.  A  concise  tutorial  on  the  application  of 
the  ANSI  procedures  for  pulsed  beam  hazard  assessment  is  given  in  Ref.  8. 

The  MPE  for  a  pulsed  beam  is  found  from  the  laser  wavelength,  the  pulse  repetition 
frequency  (PRF),  the  pulse  duration,  the  pulse  energy,  and  the  exposure  duration.  The 
laser  wavelength  and  PRF  of  the  backscattered  pulse  train  are  assumed  to  be  equal  to  the 
wavelength  and  PRF  of  the  outgoing  pulse  train.  The  term  “pulse  duration”  is  defined  in 
the  ANSI  standard  as  the  time  interval  between  the  half-power  points  on  the  leading  and 
trailing  edges  on  the  pulse.  To  find  the  width  of  the  backscattered  pulse,  one  uses  Eq.  (25) 
to  plot  the  backscattered  pulse  irradiance  at  the  center  of  the  aperture  as  a  function  of  time 
as  in  Fig.  2.  The  pulse  width  is  estimated  from  this  plot  as  the  time  at  which  the  pulse 
irradiance  is  reduced  to  half  its  maximum  value. 

The  pulsed  beam  hazard  assessment  procedure  also  uses  the  limiting  aperture  diameter, 
but  instead  of  the  effective  power  it  uses  the  effective  pulse  energy,  which  is  defined  to 
be  the  pulse  energy  transmitted  through  an  aperture  with  diameter  Df.  Exploiting  the 
rotational  symmetry  of  the  beam  model  as  in  the  derivation  of  Eq.  (32),  we  again  find  that 
the  backscattered  irradiance  from  a  single  laser  pulse  depends  only  on  the  distance  from 
the  origin.  Using  steps  similar  to  those  that  led  to  Eq.  (32),  we  find  that  the  backscattered 
irradiance  as  a  function  of  distance  from  the  origin,  p,  and  time,  r,  is  given  by 


q  ret/ 2  r+  00  r2n 

Es.p(P,t)  =  -•/  /  / 

T  Ja  JO  Jo 


l?(r,z)|2-ff(z)-Ta2(z)-Z 
[r-  cos  tp2  +  (r  •  sin  tp  —  p)2  +  z2]1-5 


■r- d< p  dr  dz. 


(34) 


A  plot  of  ESip  as  a  function  of  p  for  a  fixed  t,  e.g.,  at  t  =  t,  generated  using  Eq.  (34) 
reveals  the  spatial  properties  of  the  backscattered  irradiance  on  the  x-y  plane.  In  particular, 
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this  plot  can  be  used  to  visually  inspect  the  value,  width,  and  location  of  the  maximum 
backscattered  irradiance  on  the  x-y  plane.  If  the  maximum  backscattered  irradiance  occurs 
at  the  aperture  center,  the  backscattered  pulse  power  passing  through  the  limiting  aperture 
at  time  t,  <1 >s,/>(t),  is  found  by  integrating  Es,p(p,  /)  over  a  circular  area  of  diameter  Df 
centered  at  the  origin: 

rDf/2  _ 

<t>s,p(t)  =  2-n  /  ESip(p,  t)  ■  pdp.  (35) 

Jo 

The  effective  energy  is  then  found  by  integrating  <t>j ,p(t): 

Qs,p  =  f°°  <&s,p(t)dt.  (36) 

Jo 

A  hazardous  condition  exists  when  the  ratio  of  the  effective  energy  to  the  product  of  the 
single-pulse  MPE  and  the  area  of  the  limiting  aperture  is  greater  than  unity.  In  the  event 
that  Eq.  (34)  has  its  maximum  at  a  point  on  the  aperture  plane  other  than  the  center,  then 
the  backscattered  pulse  power  must  be  calculated  by  integrating  Es,p(p ,  0  over  a  circular 
region  of  diameter  Df  centered  at  the  point  of  maximum  irradiance.  We  note  that  the  total 
energy  in  the  backscattered  pulse  may  be  calculated  by  letting  the  upper  limit  of  the  integral 
in  Eq.  (35)  go  to  infinity  and  substituting  the  resulting  function  of  time  in  the  integral  on 
the  right  side  of  Eq.  (36). 

Further,  we  note  that  if  the  change  in  Esp(p,  t)  between  p  =  0  and  p  —  0.5  Df  is  small, 
the  backscattered  pulse  power  given  by  Eq.  (35)  can  be  approximated  by  the  product  of 
the  limiting  aperture  area  and  the  maximum  backscattered  irradiance,  Esp(0,  t).  If  we 
substitute  this  product  for  <&s,p(t)  in  Eq.  (36),  we  find  that 

n  •  Df 2  r°°  - 

Qs,p  =  —-J-  \  ES'P(0,  t)dt.  (37) 

4  Jo 

The  integral  in  Eq.  (37)  is  the  maximum  backscattered  radiant  exposure  and  is  calculated 
numerically  using  Eq.  (25).  Therefore,  when  this  approximation  is  valid,  the  ratio  of  the 
effective  energy  to  the  product  of  the  MPE  and  the  area  of  the  limiting  aperture,  which  is 
used  to  test  for  the  existence  of  a  hazardous  condition,  is  given  by  the  ratio  of  the  maximum 
backscattered  radiant  exposure  to  the  MPE 


6.  Conclusions 

High-energy  lasers  propagating  in  the  atmosphere  will  generate  backscattered  radiation 
that  may  create  ocular  and  skin  hazards  to  personnel  at  the  location  of  the  laser  source. 
The  equations  that  quantify  the  backscattered  irradiance  from  both  CW  and  pulsed  lasers 
were  derived  based  on  the  definition  of  the  backscatter  coefficient  and  beam  propagation 
theory.  The  approximate  results  derived  in  this  paper  are  appropriate  for  determining 
whether  hazardous  conditions  exist  for  personnel  in  the  vicinity  of  the  laser  source.  Our 
analysis  has  shown  that  the  backscattered  irradiance  is  a  function  of  the  laser  power  and 
the  atmospheric  backscatter  coefficient  along  the  propagation  path.  The  latter  depends  on 
the  aerosol  concentrations  and  the  atmospheric  conditions  present  at  the  test  or  deployment 
location.  Therefore,  in  order  to  assess  the  hazard  from  exposure  to  backscattered  radiation 
in  a  specific  high-energy-laser  test  or  deployment  scenario,  one  needs  to  know,  or  at  least 
be  able  to  estimate  an  upper  bound  for,  the  backscatter  coefficient  at  the  laser  wavelength 
and  the  test  or  deployment  location. 
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